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Abstract 



In a previous paper, we introduced reflection equations for interaction-round- 
a-face (IRF) models and used these to construct commuting double-row trans- 
fer matrices for solvable lattice spin models with fixed boundary conditions. 
In particular, for the Andrews-Baxter-Forrester (ABF) models, we derived 
special functional equations satisfied by the eigenvalues of the commuting 
double-row transfer matrices. Here we introduce a generalized inversion rela- 
tion method to solve these functional equations for the surface free energies. 
Although the surface free energies depend on the boundary spins we find that 
the associated surface critical exponent a s = (7 — L)/A is independent of the 
choice of boundary. 

b : 1 Introduction 

We have recently demonstrated how, by generalizing the work of Sklyanin Q , fixed 
boundary conditions may be imposed upon interaction-round-a-face (IRF) models 
whilst preserving solvability f2j. Specifically, we have shown that the Yang-Baxter 
equations and the boundary reflection equations imply commutativity of double- 
row transfer matrices. Related results have also been obtained by the authors of || 
and |§. 

Furthermore, we have found solutions to the reflection equations for the Andrews- 
Baxter- Forrester (ABF) models ||. Consideration of the fusion || |7], |^, |S| of this 
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model with fixed boundaries leads to functional relations for the double-row trans- 
fer matrices, and hence their eigenvalues. In this paper we solve these functional 
equations in the thermodynamic limit to obtain the surface free energy, from which 
we determine the surface critical exponent a s . 

In the remainder of this section we summarize the results of |2j which are needed 
for our calculation. For further details we refer the reader to this reference. In 
Section 2 we demonstrate the solution of the inversion relation for the bulk free 
energy, and then generalize this method to the surface free energy. 



1.1 The ABF models with fixed boundaries 

The ABF models || are restricted solid-on-solid models in which sites on the lattice 
take values in the set {1, 2, 3, . . . , L} subject to the condition that the values of sites 
adjacent on the lattice must differ by ±1. The Boltzmann weights depend on a 
crossing parameter A = ir/(L + 1), and a spectral parameter u. In the regimes of 
interest we have < u < A. The non-zero face weights are given by 

(a±l a \MX-u) 

\ a aT V 0i (A) 

/ a g±l\ / 1 ((a-l)A)^ 1 ((a+l)A) y /2 ^(n) 

\aTl a J V 0iOA) J t?i (A) 1 ' ' 

t a a±l\ fli(aA±rQ (1 3) 

W {a±l a J' ^(aA) " ^ 

The 0i (it) = 0i(w,p) are elliptic theta functions with nome p. We define p = exp(— e) 
with e > 0, corresponding to regime III. The critical limit is p — > 0. The boundary 
weights depend on an additional arbitrary complex parameter £, which may be 
different for the left and right boundaries. The non-zero boundary weights are 



a 



a 



0i((a± 1)A)A 1/2 M u± 0M u Ta\T0 



K a±l " = ''' ) (1.4) 



0i(aA) ) 0?(A) 



There is another form of the boundary weights which, at criticality, is independent 
of u and ^ 

K (a ± 1 ") = (^^m " 2 t^M^^A. (L5) 
V a 7 V 0i( Q A) / 01(A) 

This expression is obtained from ( |1.4| ) simply by making a complex shift in £. How- 
ever, in this paper we will take £ to be real and consider only the two forms of the 
boundary weights ( |1.4j) and ( |1 . 5| ) . From the face weights and boundary weights we 
construct a double-row transfer matrix D(w). For a lattice of width N, the entry 
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of the transfer matrix corresponding to the rows of spins a = {di, . . . , Oat+i} and 
h = {bi, . . . , &at+i} is defined diagrammatically by 

bi 61 62 63 b N bjv+l b N+1 



(a|D(«)|b> = 




The solid spins {ci, . . . , cat+i} are summed over. As the boundary weights are di- 
agonal, we must have ai = b\ and a^ + i = 6jv+i- Furthermore, these boundary 
spins, which we will call <2l and <2r, are fixed to the same values for all entries in 
the transfer matrix. The parameters £l and £r are similarly fixed for all entries. 
Defined in this way, the double-row transfer matrix exhibits the crossing symmetry 

D(A-«) = D(«). (1.6) 

More importantly, however, the double-row transfer matrices form a commuting 
family, 

D(ti)D(v) = D(v)D(it). (1.7) 

This implies that the eigenvectors of D(-u) are independent of u, so that functional 
equations satisfied by the transfer matrix are also satisfied by its eigenvalues. In 



particular, all eigenvalues of the transfer matrix satisfy the crossing symmetry ( |1.6| ). 
It should be emphasized that all the matrices in a commuting family share the same 
boundary spins and or, and the same values of £l and £r. 

The values £ = ±A/2 deserve special mention, since for these choices the isotropic 
lattice, u = A/2, has all boundary spins fixed. This is easily seen from the definition 
flOj ), as, for fixed a, only one of the choices a ± 1 gives a non-zero boundary weight. 
The non-zero boundary weights then contribute only a constant factor to each entry 
the transfer matrix. Aside from this trivial factor, the lattice exhibits pure fixed 
boundary conditions, with boundary spins alternating {a, a + 1, a, a + 1, . . . } or 
{a, a — 1, a, a — 1, ... }. 

Just as in the case of periodic boundary conditions, the face weights and bound- 
ary weights may be fused || 0, [| to form new solvable models with fixed boundary 
conditions. The functional equations which result have the same form as in the pe- 
riodic case, with the addition of some order 1 factors related to the boundary. 
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2 Functional equations 

It has been shown in that the eigenvalues of the ABF models with fixed boundary 
conditions at fusion level 1 x q satisfy the inversion identity hierarchy 



s q _ lSq+l D\u)D\u + A) = j q q ^ lS2q+1 Uf q + s 2 q Di-\u + X)D q+1 (u), 
where 1 < q < L — 1, and 

tfi(u + a) n2V 



s fc = tfi (2u+(*-l)A) 
In terms of the function 



jl = alfi, f k = (-1) 



iV 



0i (A) 



(2.1) 



(2.2) 



r-l 



on-*) = n 

J=0 



t?i(tx + (fc- j)A) 

X (A) 



(2.3) 



we can write a£ and /3£ as 



e r k (u - &)ei{u + - &)^(« + 60 



(2.4) 



= ei_ aL (u - £rX + a> + &) W« - 60^+«a(« + 60- 

If the second form of the boundary weights ( |1.5| ) is used, the $i functions in 
should be changed to $4 functions. The closure of the inversion identity hierarchy 
is governed by the conditions 



D-Hu) 



D L (u) 



0, = D L -\u) = (-l) N f- 2 a L L _ 2 /(3-2- (2-5) 



If we write the eigenvalues D q (u) with their bulk and surface terms separated 

D q {u) ~ D q h {u)D q s {u) as iV -> 00, (2.6) 

then, by virtue of the inversion relation for the fused face weights |J, D b {u) satisfies 
the functional relation 



Dl(u)Dl(u + \) = f- 1 f, 



(2.7) 



2.1 The bulk free energy 



Equation ( |2.7| ) implies that, between the inversion points u = (1 — q) A/2 and 
(3 — g) A/2, the bulk partition function per site satisfies the functional equation 



K q b (u)K q b (u + X) 



tfijX- 11)^(11 + q\) 
MW 2 ' 



(2.8) 
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This is the same inversion relation as in the case of periodic conditions, which is 
to be expected since the boundary conditions should not affect the bulk behaviour. 



Equation (|2.8|) has been solved previously |M 0, [LOj , but we include the solution here 



for completeness, and as an introduction to the generalized methods that follow. We 



use the standard techniques developed by Baxter PHI . The assumption that the so- 
lutions are analytic between the inversion points, and that they may be analytically 
continued a small distance outside the strip, along with the relation (|2.8|) and the 
crossing symmetry 

D q {{2 - q)\ - u) = D q {u) (2.9) 

uniquely determines the free energies. The assumption of analyticity may be justified 
by studying the zeros of the largest eigenvalue D q (u) for large finite N. In the critical 
case, it is seen that for the fusion level 1 x q, the strip —q\/2 < Re(w) < (4 — q)X/2 is 
free of order N zeros JTJJ. Furthermore, there can be no poles inside this strip, since 



for finite N the independence of the eigenvectors on u implies that the eigenvalues 
are simply linear combinations of products of Boltzmann weights. Since none of 
the Boltzmann weights have poles, neither can the eigenvalues. This assumption of 
analyticity implies that the logarithms of the partition functions may be expanded 
in a Laurent series in powers of exp(27ru/e), 

oo 

In K q h (u)= J2 c k e 2k ™ /£ . (2.10) 

k=— oo 

We rewrite the right hand side of (|2.8|) using the "conjugate modulus" transforma- 
tion, which, in terms of the function 

oo 

E(x,p) = JJ(1 -p^x)^ -p n x~ l )(l -p n ), (2.11) 



n=l 



is given by 

tfi(u,p)=(jj e -{u-*m 2 /s E (e' 27Tu/£ , q 2 ) (2.12) 

4 ( u ,p) = Qy /2 e -(W2) a /^(_ e -W« >§ 2) (2.13) 

where p = exp(— e) and q = exp(— 7r 2 /e) are conjugate nomes. With both sides 
of Q2.8|) expanded in powers of exp(27ra/e), we match coefficients and impose the 
crossing symmetry (|2.9|) to obtain the solution 



o, x / x cosh[(7r - 2\)irk/e 



k- 
oo 



k=l 



J ksmh.(ir 2 k/e) 

k smh.(n 2 k/e) cosh(A7r/c/e) 



Ecosh[(7r — (q + l)A)7rfc/e] cosh[((2 — q)X — 2u)nk/e] , ^ 
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where 



co(u) = ^ [(tt - gA)(g - 1)A + 2u((2 - g)A - u)] . (2.15) 



The analytic continuation of this function to the region —q\/2 < Re(tt) < (4 — q)X/2 
gives the bulk behaviour of the largest eigenvalue of the transfer matrix 

D l( u ) ~ [«b(«)] 2iV as iV ^ oo. (2.16) 
From ( p.14 ) it is easy to show that inside the interval —q\/2 < u < (2 — q) A/2, 

K-fV + AK + V) = + A) wheng>l. (2.17) 

When g = 1, it can also be shown that inside —A/2 < u < A/2, 

^ _1 (u + A)^ +1 (w) tt ^ sinh(A-2|M|)7rfc/e 

^ --(A - 2M) - ^ fccosh(A7rfc/£) ' (2- 18 ) 



In 

so that 



K q h (u)K q h (u + X) 



|<~V + AK + V)| < |«g(u)«£(u + A)| when g= 1. (2.19) 
Putting together equations fl2.8Q , ([2.16D , ( |2.17| ) and ( |2. 19| ) , we therefore obtain 

nm ( Dl 1 (u + X)Dl +1 (u) \ = |0 when q = 1, 
f-i(u)f q (u) ) \l wheng>l, 

which is consistent with the critical bulk behaviour described in [JE^]. We need this 
result for the derivation that follows. 

2.2 The surface free energy 

Recalling the separation of the bulk and surface terms (|2.6|) , and using the functional 
equation ( |2.7| ) and the bulk behaviour of the eigenvalues ( |2.20| ), the inversion identity 
hierarchy in the thermodynamic limit becomes, for —qX/2 < Re(tt) < (2 — q)X/2, 



s q „ lSq+1 DI(u)D«(u + X) 



7oS-iS 3 if q = 1, 

7*_ 1 s_ 1 s 2(?+1 + s 2 q DI-\u + X)DI +1 (u) otherwise. 

(2.21) 



In the case of the unfused model, q = 1, we therefore have an inversion relation 
which allows us to calculate the surface free energy^. Explicitly, 

, , , tf 1 (2A-2«)tfi(2A + 2u) r. , R/ , 



3 In the general case of p x q fusion, an analogous functional equation is derived when q = p. 
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where e L (w) and e R (w) are defined by, for a and £ corresponding to the left or right 
boundary as appropriate, 



In the unfused case, the crossing symmetry of the eigenvalues implies that 



The solution of the inversion relation for the surface free energy proceeds in a similar 
fashion as that for the bulk free energy, although we must justify the assumption 
of analyticity separately. There can of course be no poles between the inversion 
points for the same reason as there are none in the bulk. However, we are now con- 
cerned about order 1 zeros inside < Re(it) < A. Certainly zeros occur on the line 
Re(u) = A/2, but, guided by the derivation of conformal weights for periodic bound- 
ary conditions ||12| , we associate these zeros with finite-size (order 1/N) corrections 
rather than surface effects. 

Numerical studies show that zeros do occur on the real u axis inside the strip 
— A/2 < Re(tt) < 3A/2, but never inside the strip < Re (it) < A. For certain values 
of £, zeros occur at the inversion points u = and A. These values are determined 
by the zeros of e L,R (w), and are found to be £ = 0, —aX and (L + 1 — a)X. For all 
other values of f (up to periodicity), the interval between the inversion points is free 
of zeros. We therefore conclude that the surface partition function per site, k s (u), 
is analytic in this region. 

With this assumption of analyticity, and the imposition of the crossing symmetry, 
the solution of ( p. 22] ) is 



e L - R (n) = el(u- oel(u + ooUu - i)e l a (u + §. 



(2.23) 



Ks(u). 



(2.24) 




sinh[(7r — 3\)nk/e] sinh(A7rA;/£) cosh[2(A — 2u)7rk/s] 




(2.25) 



+ (tt - 3\)\/e + In k^(u) + In nf(u), 
where In and In k r are given for generic a and £ by 




CO 



cosh[(aA + £ — |£|)7rA;/e] cosh[(7r — a\ — £ — |£|)7rfe/e] cosh[(A — 2u)nk/e] 



2 E 




k=l 



and 



c (it) = - [(aX 20 + (|f | - 2A)tt + (2 - a 2 )A 2 + 2u(X - «)] . (2.27) 
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In deriving this expression, £ is assumed to satisfy the inequality 



-7T < (1 -a)X < £ < (L-a)X < 7T. (2.28) 

We note that, as one would expect, the height reversal transformation a — > L+ 1 — a 
and £ — > — £ leaves ( |2.26| ) unchanged. If the second form of the boundary weights 
(P-.5Q is used, each term in the sums of ( |2.26|) should be multiplied by (— l) k , which 
alters the critical behaviour. 

The temperature variable t is identified in || to be given by t = p 2 . The 
behaviour in the critical limit t — > + is found by applying the Poisson summation 
formula to the above expression for \nn s ( 2.25[) . We find that the leading-order 
singularities of In k s have the form 

f^~K> lfL "°° rl (m ° d4) ' ^9) 
J S \W 4A lnt if Lee 3 (mod 4). V ; 

When L = 2 (mod 4), In k s is regular. The leading singularity of In k^ ,r is in general 
of higher order than that of lnK s . In addition, In Kg' is regular in the following 
situations: 



£ = k\, k e Z 
L even ^ £ > and a odd L odd 

£ < and a even 



f = (2fc + l)A/2, fceZ 
a odd 



These exceptions aside, In k l ' r has the leading-order singularities 

f L,n U n/2X if Lis even, 

/sing 1 W 2A lnt if Lis odd. 1 ' ] 



When the second form of the boundary weights (|1.5|) is used, the leading-order 
singularities are, aside once again from the above exceptions, 

,l,r JV /4A if L is even, 

/stag ~ | r /4A lnt if L is odd. (2 ' 31) 

In this case the function In k^ ,r vanishes at criticality. 

The surface critical exponent a s is defined in analogy with the bulk critical 



exponent a |13J, |14 



with the specific heats 



C h ~\tf a C a ~\t\- as , (2.32) 
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Our result for the surface critical exponent a s is therefore 

a s = 7 -^. (2.34) 

The ABF model with L = 3 corresponds to the two-dimensional Ising model. In 
this case we have 

/sing -tint, (2.35) 



in agreement with the calculations of McCoy and Wu |15] for free boundary con- 
ditions, but with a lattice rotated by 45° with respect to the one considered by 
them. 



3 Conclusion 

From the inversion identity hierarchy, and from the known solution of the bulk 
free energy, we have derived an inversion relation for the surface free energy of the 
ABF models with fixed boundary conditions. We have solved this inversion relation, 
subject to justifiable analyticity assumptions, and thus obtained an expression for 
the surface free energy. Finally, we have analysed the critical behaviour of the 
surface free energy to obtain the surface critical exponent a s . 

In this paper we have considered the bulk form of the free energy and its surface 
(order 1) correction. In a future publication we will study the finite-size (order 1/N) 
corrections at criticality, and hence derive the central charges and conformal weights 
of the ABF models with fixed boundaries. 



After this work was completed, we received the preprints |16| and [17| , in which 
the authors derive surface critical exponents of the eight-vertex and ABF models 
using methods similar to ours. 
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